Polymers with attractive interactions on the Husimi lattice 



O 
O 
(N 

1—5 



Pablo SerrfQ 

Facultad de Matemdtica, Astronomia y Fisica 

Universidad Nacional de Cordoba 

Ciudad Universitaria - 5000 Cordoba 

Argentina 

Jiirgen F. StilclQ 

Instituto de Fisica 

Universidade Federal Fluminense 

Av. Litordnea s/n 

24210-340 - Niteroi, RJ 

Brazil 






-)— > 

a 

o 
o 



> 

m 

\o 
en 

o 

o 

-I— • 



I 

O 

o 



X 



Welchy L. Cavalcanti and Kleber D. Machadcp 

Departamento de Fisica 

Universidade Federal de Santa Catarina 

88.040-900 - Florianopolis - SC 

Brazil 

(Dated: February 2, 2008) 

We obtain the solution of models of self-avoiding walks with attractive interactions on Husimi 
lattices built with squares. Two attractive interactions are considered; between monomers on first- 
neighbor sites and not consecutive along a walk and between bonds located on opposite edges of 
elementary squares. For coordination numbers q > A, two phases, one polymerized the other non- 
polymerized, are present in the phase diagram. For small values of the attractive interaction the 
transition between those phases is continuous, but for higher values a first-order transition is found. 
Both regimes are separated by a tricritical point. For q=4 a richer phase diagram is found, with an 
additional (dense) polymerized phase, which is stable for for sufficiently strong interactions between 
bonds. The phase diagram of the model in the three-dimensional parameter space displays surfaces 
of continuous and discontinuous phase transitions and lines of tricritical points, critical endpoints 
and triple points. 

PACS numbers: 05.50.-|-q, 61.41. -|-e, 64.60.Ht 



INTRODUCTION AND DEFINITION OF 
THE MODEL 



Self-avoiding walks have been found to be useful mod- 
els for the study of the behavior of polymers for quite 
a long time [l|. The self- avoidance constraint in general 
makes these models difhcult to solve on regular lattices. 
On lattices with hierarchical treelike structure, however, 
like the Bethe |3| and the Husimi [3 lattices, it is not 
difficult to solve such models exactly. From the point of 
view of critical phenomena, these solutions lead to clas- 
sical or ideal chain critical exponents, but non-universal 
features of the phase diagram may be closer to the ones 
observed on regular lattices than those provided by the 
usual mean- field methods Q. 

The effect of monomer-solvent interactions on the be- 
havior of polymers diluted in poor solvents may be in- 
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eluded in the model by allowing attractive interactions 
between segments of the chains. This induces a com- 
petition between repulsive, excluded volume interactions 
and the attractive short range interactions. As these 
latter interactions become sufficiently strong, the chains 
may change from an extended to a collapsed state |a]- 
While this collapse transition, usually identified with the 
0-point, appears as a tricritical point in mean-field ap- 
proximations and non-classical approximations on three- 
dimensional lattices, in two dimensions the situation does 
not seem to be so simple. Transfer-matrix calculations [g 
and exact Bethe-ansatz results Q for a 0{n) model with 
four-spin interactions on the square lattice lead two phase 
diagrams where the second order transition line between 
the polymerized and non-polymerized phases ends at a 
multicritical point whose precise nature is not clear from 
these calculations, but which is definitely not a tricriti- 
cal point. In the limit n — > this model corresponds to 
self-avoiding walks with attractive interactions between 
bonds of the walk which are located on opposite edges 
of elementary squares of the lattice. The four-spin inter- 
actions in the magnetic model are related to interactions 
between bonds on the corresponding polymer model. On 
the other hand, studies of the behavior of self-avoiding 



walks on the square lattice with attractive interactions 
between monomers located on first-neighbor sites but not 
consecutive along the walk point to a tricritical collapse 
transition 'y| . The main motivation of this work was to 
address the question of apparently quite similar models of 
attractive walks (bonds or monomers interacting) leading 
to qualitatively distinct phase diagrams. The problem 
of models with interactions between monomers or bonds 
only on Husimi lattices was considered in the literature 
0,113) [III- On a four coordinated (g — 4) Husimi tree, 
the phase diagram of the model with monomer- monomer 
interactions is qualitatively similar to the one found in 
general when g > 4. In the parameter space defined by 
the activity of a monomer (x) and the Boltzmann factor 
of the elementary interaction between bonds (w), a non- 
polymerized phase is stable at low values of x, whereas 
a polymerized phase is found at higher activities. The 
transition between those phases is continuous at low val- 
ues of Lu, but becomes discontinuous as lu is increased. 
These two regimes are separated by a tricritical point. 
When the interactions are between bonds (Boltzmann 
factor k) , a third phase is stable in part of the parameter 
space. This phase is a dense phase (all sites are visited 
by the polymer) and the transition continuous transition 
line between the non-polymerized and the polymerized 
phases ends at a critical endpoint. The phase transi- 
tion between the polymerized and the dense phases may 
be continuous or not, a tricritical point being found on 
this transition line. Transfer matrix calculations for this 
model on the square lattice suggest that this picture may 
be observed there as well ^^■ 

We solve a model of interacting self-avoiding walks on 
Husimi lattices (core of the Husimi trees), built with 
squares Q . At each site of the lattice the ramification 
of squares is equal to a, and therefore the coordination 
number will be q = 2{a + 1). Solutions of models on 
such lattices may be considered approximations to the 
solution on hipercubic lattices of the same coordination 
number, so that a Husimi lattice with a = 1 leads to 
a solution which approximates the one on a square lat- 
tice, the solution for a = 2 may be an approximation 
for the solution on the cubic lattice, and so on. The 
Husimi lattice solution may be considered to be the third 
member of a sequence of approximations whose first two 
are regular mean-field and Bethe lattice solutions. In a 
mean-field calculation no correlations are taken into ac- 
count, while Bethe lattice and Husimi lattice solutions 
short range correlations are taken care of. It should be 
stressed, however, that all these solutions lead to classi- 
cal critical exponents. In some cases, it has been shown 
that Bethe and Husimi lattice solutions may show fea- 
tures of the phase diagram of models on regular lattices 
which are not present in the corresponding mean-field 
approximations 0,EQE3- 

The model we consider is of self-avoiding walks on a 
Husimi tree, with the initial and final monomer of each 
walk located at the surface of the tree, so that the density 
of endpoint monomers in the core of the tree is always 
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FIG. 1: A Husimi tree with a — 1 and three generations with 
two polymers on it. The statistical weight of this configura- 
tion is equal to x^^co^k^. 



equal to zero. We associate an activity x to each bond of 
a walk, and include an interaction energy em for each pair 
of monomers on first neighbor sites of the lattice with no 
bond of the walk between them. Also, an interaction 
energy £(, is associated to each pair of bonds located on 
opposite edges of elementary squares. The grand parti- 
tion function of the model on a lattice with N sites may 
be written as 



Y{X,LU,K-N)^Y1 



X^^LO^- 



.Nil 



(1) 



where the sum is over all configurations of the walks on 
the lattice, Ni, is the number of bonds in the configura- 
tion, Nim is the number of pairs of interacting monomers, 
and Nil) is the number of interacting bonds. The Boltz- 
mann factors which correspond to the interactions are 
given by w = exp(— em/^BT) and k = exp(— et/ZcsT). A 
configuration of walks and the corresponding statistical 
weight may be found in Fig. ^ 



II. SOLUTION OF THE MODEL ON THE 
HUSIMI TREE 

To solve the model on the Husimi lattice we use a re- 
cursive procedure, defining subtrees of the Husimi tree 
and establishing recursion relations between the partial 
partition functions of the model on the subtrees, for fixed 
configurations of the root site. Fig. |21 shows the three 
possible root configurations of a subtree, labeled by the 
number of bonds incident at the root site from above, 
and a diagram illustrating how to obtain the partial par- 
tition functions of a (n-f l)-generations subtree from the 
partial partition functions of 7i-generations subtrees. Ini- 
tially, we consider three subcases for the root configu- 
ration with no incident bond, defined by the number of 
monomers in the first neighbor sites to the root site (0, 
1, or 2). To obtain the recursion relations for go.Oj 5o,ij 
<7o.2, 5i, and ^2, we consider all the possibilities of attach- 
ing three sets of a n-generations subtrees to the vertices 



of the elementary square at the root of the new {n + 1)- 
generations subtree. The recursion relations are obtained 
so that the activity of the bonds and the Boltzmann fac- 
tors of the interactions between monomers and bonds in 
the elementary square at the root of the new (n + 1)- 
generations subtree are considered in the iteration. We 
then notice that the partial partition functions for the 
configuration with no incident bond appear in only two 
combinations in the recursion relations, which are: 
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5o = 50.0 + 50, 1 +5o,2 



(2) 



and 



(a) 



53 = 50,0 + ^.90,1 + ^^50,2- 



(3) 



If we call gi, i = 0, 1, 2, 3 the partial partition functions 
of the model defined on an n-generations subtree and g^ 
the same functions on a n + 1-generations subtree, we 
may write the recursion relations as 



5o 
5i 

9's 



where 



and 



SHg^'^ + (1 + 2uj)H^goa 



2xF^g^ + Lu^H^ + 2xu;F^H + x'^F^g^, 
2xF{gl'' + 2ujHg^ + u;^H^ + xg^g^ + 
xuj'^Hgl + x^ungl" + xlu'^kF'^), 
x^F^{g^ + uj^H + 2xu;Kg^), 
.g^-^ + Sffgo^- + 2u;{Hgl'^ + 
LoH^g^ + xF^g^)+uj^{H^g^ 
2xujF^H + x^F^g^; 



uj'H^ 



F = cr5i5^ 



(<T-1) 



H = cr.g2.g3 



(a-l) 



a{a 



1) 2 . 
— 5i53 



(^-2) 



and c : 



53 
50' 



It is convenient to define the ratios 

51 , .92 
a = — , &= — , 

50 50 

which obey the following recursion relations 

a' = 2x/(l + 2ujh + uj^h^ + xc" + 

xu^c^h + x^ujuf?" + xup-K,f)lq, 
b' = x^f{l+uj'^h + 2xujKc'')/q, 
c' = {l + h + 2ujh + 2,Lj'^h^ + 2ujx,f + 
uj^h^ + 2xu;^fh + x^uj^c" f)/q, 



where 
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1 + 3/i + (1 + 2Lj)h^ + 2xf + 



2n f2 



(4a) 

(4b) 
(4c) 



(4d) 



(5) 



(6) 



(7) 



uj^k' + 2xuj}^h + x^c''f 
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FIG. 2: (a) The possible configurations of the root site of a 
subtree, (b) A 3-generations subtree is built attaching three 
2-generations subtrees to the vertices of the new root square. 
In this example a = 1. Two 2-generation subtrees have root 
configuration 1, while the other one has root configuration 0,0. 
The resulting root configuration of the 3-generations subtree 
is 2. 



The partition function of the model on the Husimi tree 
may then be obtained if we consider the operation of at- 
taching four sets of a subtrees to the central square of the 
tree. In Fig. |3|the contributing configurations of the cen- 
tral square are shown in the order of the corresponding 
monomials appearing in the resulting expression below 



y„(a;,w,K) 



^4cr 



5d" + H^H + 2gt,''H' + 
^xgl'^F^ + Auj^g^H^ + 
Sxujg^F^H + ix^g^g^F^ 



Acog'o^H 



iu^H' 



Axoj^F^H^ 



2kx^uj^F'^ 



Ax^uj^g^F^H + Ax^ujngfF^. 



(10) 



We expect the thermodynamic behavior of the model on 
the Husimi tree to be quite different from the one found 
on regular lattices, since the surface sites dominate in 
the thermodynamic limit, when the number of iterations 
n — > 00 [lq |. We therefore will focus our attention on 
the behavior in the central region of the tree, which we 
will refer to as the Husimi lattice J3|. Considering the 



contributions to the partition function in Eq. ^J we 
may calculate the mean numbers of bonds, monomer- 
monomer interactions and bond-bond interactions in the 
central square of the tree, which are given by 



Pb = 4xf{l + 2ujh + 2xc'' + 

uj^h^ + Kxuj'^f + 2xu?hc'' + 

ix^UJKC^'')/d, 

tmm = 4:Uj{h'^ + 2ujh^ + 2xph + 

LU^h^ + ixU^fh^ + KX^Ujf + 

2x'^ujfhc'' + x^Kpc^")ld, 
Pbb = 2x^ujKf{Lj^f + 2xc^'')/d; 



(11a) 



(lib) 
(lie) 



where 



d = l + Ah + 2h'^ + Aioh^ + Axf + Auh^ + 
^xujfh + Ax'^c^2'^ + uj%^ + Axuj^fh^ + 
2KX^Lj^f^ + 4:X^uj'^c''fh + Ax^ujnc^^f, (12) 

calculated at the fixed point of the recursion relations 
(Eqs.EJ. 

The thermodynamic behavior of the model is deter- 
mined by the fixed points of the recursion relations Eqs. 
IHl each of which correspond to a thermodynamic phase. 
We investigated the stability regions for each of the 
fixed points wc found, which are three, in general: (a) 
a = 6 = 0, c = 1, the non-polymerized fixed point, cor- 
responding to pb = Pmm = Pbb = 0, (b) a, 6, c 7^ and 
finite, which is the regular polymerized fixed point, where 
nonzero densities are found, and (c)a —> c», and &, c ^ 
and finite, which we call saturated polymerized phase, 
since pb — Pmm ~ 2, pbb = 1 in this phase. This latter 
is stable in a region of the phase diagram only for the 
four-coordinated Husimi lattice a = 1, being absent in 
the phase diagrams of higher coordinated Husimi lattices. 
To study the stability region of the saturated phase it is 
convenient to rewrite the recursion relations [S] in terms 
of the following new variables: 



« 90 

a = —, 

91 

/3 = ^, 
91 

93 

91 



and 7 



(13a) 
(13b) 
(13c) 



In these variables, the saturated fixed point is located at 
the origin (a = /? = 7 = 0). 

The stability regions of phases (a) and (c) may be 
found analytically, due to their simplicity. At the stabil- 
ity limit of the non-polymerized phase the largest eigen- 
value of the jacobian associated to the recursion relations 
IHl calculated at the fixed point (a, 6 = and c = 1) is 
equal to unity. The result is: 



K < 



1 — 2xa — 2x a 



while the stability region of the dense polymerized fixed 
point is obtained in a similar way using the recursion 
relations written in the other set of variables El One 
obtains {a ~ 1): 



X > 



-1 + 8k + Lu - 8kuj + Ak^uj 



(15) 



It should be remarked that for n = \ + \/3/4 f« 1.86 
the second member of the inequality above diverges, so 
that the dense phase is never stable for k < 1 -f ^3/4. 
In particular, this is true for the model with interactions 
between monomers only {k = 1), as remarked by Pretti 
[lOl . The stability limit of the regular polymerized phase 
(b) may be found numerically. 

The critical surfaces in the phase diagrams will eventu- 
ally end at tricritical lines. These lines may be obtained 
requiring the corresponding solution to be a double root 
of the fixed point equations. Some algebra furnishes the 
tricritical condition for the (a)-(b) critical surface: 



P{uj,x) 



= —uj + 7x — 2ujx — 16a; -(- 
A{uj - 2)uj x^ + 8{l + {uj 



lOujx-' - 

- 1)lu)x'^ + 



2{1 + 2lj{u; - l))x^ 



0. 



(16) 



2x^auj 



(14) 



A similar calculation for the tricritical condition on the 
critical surface where phases (b) and (c) are equal leads 
do: 

r{uj,K) = -2{uj-lf + {-l+uJf{63 + 2uJ + u^)K- 
S{99 - ISSw -I- 85w^ -f 4cj'*)k2 + 
16(311 - 553cj + 2yiu? -I- 25w'^)k^ - 
32(498 - 786w + 213^^ -|- 76w'')k'' + 
16(1483 - 1876w - 44^^ + 454w'')k5 - 
128(112 - 120cj - 65cj2 + 1<6u:^)k^ + 
128(9 - 22uj - 42w^ + 50cj'*)k'^ - 
\Q2\u?{-\ + 2<J)k^ + 
256u;''k^ = 0. (17) 

To obtain the phase diagrams of the model, in the pa- 
rameter space defined by x, w, and k, we find which fixed 
point is stable at each point of the parameter space. Sur- 
faces of this space where the stability limits of two phases 
are coincident are critical surfaces, and regions where 
more than one fixed point is stable are related to first 
order transitions. To obtain the location of these tran- 
sitions one may use a Maxwell construction, although 
it is sometimes possible to find the free energy of the 
model on a treelike lattice using appropriate recursion re- 
lations 4] . Due to the simplicity of the non-polymerized 
and saturated phases, the partition function per elemen- 
tary square of the lattice may easily be calculated, and 
this result simplifies considerably the determination of 
the coexistence surface of these two phases. Since in the 
non-polymerized phase only configuration 1 of figure 13 is 
present in the core of the lattice, the partition function 
per elementary square will be y„p = 1. In the saturated 
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FIG. 3: Possible configurations of the central square of the tree. Each of them contributes with a monomial in the calculation 
of the partition function of the model on the Husimi tree. 



phase, configuration 11 will dominate, and therefore for 
this phase we have Zg = Inx^uP'. Thus, for the four- 
coordinated lattice the coexistence surface, in the region 
where phases (a) and (c) are stable, is given by: 



2kx^u? = 1. 



(18) 



Since the phases (a) and (c) have different densities, the 
transition between them is always of first order, as may 
be verified in the phase diagrams below. 



III. PHASE DIAGRAMS 

As mentioned before, for cr > 2 only the non- 
polymerized and the regular polymerized fixed points of 
the recursion relations (Eqs. IHl are stable, and there- 
fore no dense polymerized phase appears in the phase 
diagram. The non polymerized phase is stable for small 
values of the activity x, and as the activity is increased, 
eventually the regular polymerized phase becomes sta- 
ble. The stability limits of both phases are coincident at 
low values of the interaction (lo and k close to one), and 
thus the phase transition between them is continuous. As 
the strength of the interaction is increased, however, the 
transition will become discontinuous, and thus the crit- 
ical surface (where the both phases have the same den- 
sities) is separated by a tricritical line from the surface 
where both phases coexist. This phase diagram, where 
the collapse transition corresponds to a tricritical point, 
is expected for this problem since the theta point was 
recognized as a tricritical point in the pioneering work 
by de Gennes (J3])j and thus we will concentrate now on 
the case of the four coordinated lattice. 

For cr = 1 the dense phase is stable in part of the 
phase diagram, and therefore richer phase diagrams are 
obtained. We will consider here constant u) cuts of the 
phase diagram. Perhaps the most interesting constant k 
diagram is the one with interactions between monomers 
only [k = 1), which may be found in figure 2 in the com- 
ment by Pretti [i3. Essentially three different types of 
phase diagrams are found in the (x, k) plane for increas- 
ing values of u. 
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FIG. 4: Phase diagrams of the model for a = \ and u) = \ 
Continuous lines are second order transitions and dashed lines 
are first order transitions. Tricritical points are indicated by 
circles and critical endpoints are represented by squares. 



1. For 1 < oj < wi the critical polymerization line 
ends at a critical endpoint located at the conflu- 
ence of the coexistence lines of the dense phase with 
the other two. The particular case with interac- 
tions between bonds only (lo — 1) is depicted in 
figure 2lThe dense phase and the regular polymer- 
ized phase are separated by a transition line which 
may be of first or second order, a tricritical point 
separating these two cases. The critical endpoint 
becomes a tricritical point at iO\ ~ 1.15301, k\ ~ 
4.46985, xi ~ 0.29007. These values may be ob- 
tained noting that this point in the parameter space 
is located on the coexistence surface of the non- 
polymerized and the dense phases (equation ^J), 
and on the tricritical line, defined by the stability 
limit of the non-polymerized phase ('eauation ll4l as 
an equality) and the tricritical condition (equation 
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FIG. 5: Phase diagrams of the model for cr = 1 and tj = 1.18. 
The triple point is indicated by a triangle. 



2. For LJi < Lo < UJ2 tricritical points are present on 
the boundaries of the polymerized phase with both 
the non-polymerized and the dense phases. Three 
coexistence lines meet at a triple point, as may be 
seen in figure [S] As w is increased, the tricriti- 
cal point on the boundary between the polymer- 
ized and the dense phases moves closer to the triple 
point. These two points meet, becoming a critical 
endpoint, at 0J2 ^ 1.21717, K2 ^ 4.08985, X2 ^ 
0.28726. 

3. For Lo > LU2 the transition between the polymer- 
ized and dense phases is always continuous, and 
this critical line ends at a critical endpoint. An ex- 
ample is shown in figure [^ The boundary between 
the non-polymerized and the polymerized phases 
displays a tricritical point if w < 1.54508 and is 
always of first order for higher values of lo. 

For the determination of the first order boundaries in- 
volving the regular polymerized phase, a Maxwell con- 
struction was done using the pair of conjugated variables 
K and pbb- Other options are possible, this one was cho- 
sen for simplicity. The result of this calculation is not 
expected to depend on the choice of the variables. For a 
similar model the Maxwell relations were tested explic- 
itly [131 and also the free energy was obtained directly 
using a proper iterative procedure |j]. 



IV. CONCLUSION 

The solution of models for polymers with attractive 
interactions on Husimi lattices built with squares leads 
to the expected phase diagrams when the ramification 
of the lattice a is equal or larger to 2. In this case, 
only one polymerized phase is found, separated from the 
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FIG. 6: Phase diagrams of the model for a — 1 and 
uj — 2. Since uj > 1.54508 the transition between the non- 
polymerized and the polymerized phase is always of first or- 
der. 



non-polymerized phase by a first- or second order tran- 
sition line. The two lines are separated by a tricritical 
point, which is associated to the collapse transition of the 
polymers. At a four-coordinated Husimi lattice {a = 1), 
however, a second polymerized phase is stable at high 
values of the Boltzmann factor k of the interactions be- 
tween bonds. Since in this phase all sites of the lattice 
are incorporated into the polymers, we called it dense 
polymerized phase. 

The phase diagrams obtained for the model on the 
four coordinated Husimi lattice may offer an explana- 
tion for apparently confiicting results in the literature re- 
lated to the collapse transition of polymers on the square 
lattice. Transfer matrix and finite size scaling calcula- 
tions of a model with interactions between monomers on 
first neighbor sites 8] lead to compelling evidences that, 
as UJ is increased, the critical polymerization line ends 
at a tricritical point, but exact Bethe ansatz arguments 
for a magnetic model which is equivalent to a polymer 
model with attractive interactions between bonds do not 
indicate a tricritical collapse transition point Q- In the 
Husimi lattice solution of the model presented here, the 
collapse transition point is a tricritical point for the case 
of interacting monomers and a critical endpoint when 
the interaction is between bonds. This is consistent with 
the results known for these models on the square lattice. 
It is possible to obtain many features of the problem of 
directed polymers with interacting bonds on the square 
lattice exactly [3|- Although the collapse transition for 
this case is also a multicritical point, the details of the 
phase diagram are quite different from the ones found 
here. A rather unphysical feature of this model is that a 
polymerized phase of zero density of monomers occupies 
a finite region of the phase diagram, even in the absence 
of any attractive interactions. 



It is of some interest to find out how the quahta- 
tively different phase diagrams for interactions between 
monomers and bonds change into each other. Thus, we 
considered here a more general model where both inter- 
actions are present. We then found that the transition 
between phase diagrams were the collapse transition is 
a tricritical point (interactions between monomers) and 
those were it is a critical endpoint (interaction between 
bonds) occurs through an intermediate phase diagram 
were two tricritical points are present and three coexis- 
tence lines meet at a triple point. It should be stressed 
that the two tricritical points found in the phase diagram 
are not part of the same tricritical line in the full pa- 
rameter space (a;, w, k). This poses the question if those 
two distinct tricritical points actually exist for the model 
defined on two-dimensional lattices and, if this is true, 
if they share the same set of tricritical exponents. Al- 
though experimental studies of the collapse transition of 



polymers confined in a two-dimensional surface have been 
done ,19.] and seem to confirm the tricritical nature of the 
transition, monodisperse polymer solutions were used for 
them, while in our calculations the chains are polydis- 
perse. 
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